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ABSTRACT 

We investigate the sensitivity of the semileptonic processes e+e" ^ i'ui q (f, i = 
e or ^, on the non-standard trilinear gauge couplings, using the optimal observables 
method at Linear Collider energies. Our study is based on the four-fermion generator 
ERATO. Taking into account all possible correlations between the different trilinear 
gauge coupling parameters, we show that they can be measured with an accuracy of 
10"'^ to 10~^ for typical Linear Collider energies and luminosities. 
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The future e~^e~ linear colliders (LC), with energies ranging from a few hundreds of GeV up 
to a couple of TeV, provide particle physics with an enormous potential for studying Nature to 
the deepest level ever achieved fll. Although at these energies direct searches for new particles 
and their interactions will eventually attract most of the physics interest, LC offer also the unique 
possibility to study to an extremely high accuracy, the properties of the existing particles, like 
those of the massive electroweak gauge bosons W and Z. Therefore, an important project at the 
LC energies will be the determination of the trilinear gauge couplings (TGC) |2|, ^, which are 
a characteristic manifestation of the underlying non-Abelian symmetry of elementary particle 
interactions Q and at the same time an interesting probe of New Physics (NP). 

In order to study the TGC we need a parameterization of the vector gauge boson interactions 
that goes beyond the Standard Model . The most general such parameterization is given by 

Ctgc = E -i(^9vww{9Y{Vf,W-''''W+ -V^W-^>'^'W-) + KvV^,W+>'W- 
v=■y,z 
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where 



and 



In Eq.(ll) is the M^-boson field, and the usual definitions g^ww = Ij dzww = ctg Oy^ are 
used. In the Standard Model we have gj = gf = I, k-^ = kz = while all the other parameters 
are vanishing at tree level. In searching for possible TGC, it is more convenient to express them 
in terms of their deviations from the Standard Model values. For this we define the deviation 
parameters [||, ^: 



= {gi - l)ctg 6w Xj = Kj -1 xz = {kz - l)ctg By, - 6z , (2) 

while we throughout assume gj = 1, disregarding the possibility of an anomalous contribution 
to the electromagnetic form factor of W^. We note that the NP contribution described by the 
interaction Lagrangian in Eq.(|lj), becomes linear when expressed in terms of the above deviation 
parameters and A^, Xz, as well as the C- and P-violating couplings. 

During the last years, considerable progress has been achieved concerning the understanding 
of the physics underlying the non-standard boson self-couplings. As showed in reference |^, 
the deviations from the Standard Model TGC couplings can be parameterized in a manifestly 
gauge-invariant way by using the effective Lagrangian approach and considering gauge-invariant 
operators involving higher-dimensional interactions among the gauge bosons and the Higgs field. 
These operators are scaled by an unknown parameter A^vp describing the characteristic scale of 
some high energy New Physics, generating at low energies the effective interaction Ctgc as a 
residual effect. In order to generate all kinds of TGC appearing in Eq.(|^), we need operators with 
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dimension up to twelve. On the other hand, restricting ourselves to SU{2)i x ?7(l)y-invariant 
operators with dimension six, which are the lowest order ones in a l/A^vp expansion, we end up 
with the following list of operators capable to induce TGC NP couplings |^ |9|, 11, 14 1: 



and[| 



Ow = ^{W^p^WP,)-W'^^ (3) 



Obw = ^^'^■W^'^B^, 



Ow = l^_{W^,xW^,)-W\, (4) 



where 

1 l,l,nrr „ TTT-A"^ ^ 



Bt^- = -e>^-P-Bp^ , W'^=-e>''""'Wp^. (5) 
In Eq.(^) and Eq.(§), Tj describe the Pauli matrices, 

B^u = d^By - dyB^ 

is the U{1)y gauge field strength, 

w^, = d^w, - d,W^ - gW^ X 
is the field strength for the SU{2)l gauge field W^, and the Higgs doublet is written as 



\ V2 



1 („ + i7 + i<^0) 



while Dfj^ is the covariant derivative 



D^ = df, + ig--W^ + ig'YB^ , 

and Y is the hypercharge of the field on which Z)^ acts. Finally e = gsin^^ = g' cos Ow 

In the list of Eqs.(|^,^), we have included all dini=6 purely bosonic operators contributing to 
the trilinear gauge interactions, except those which give also a tree level contribition to LEPl 
observables, (and of course those which give no TGC at all). This constitutes part of a consistent 
general strategy for searching for any purely bosonic NP interaction. According to this strategy, 
the measurement of TGC provides the most efficient way to study the operators appearing in 
Eqs.(^,^, while the rest of dim=6 operators can be most efficiently disentagled and constrained 
either by high precision measurements (LEPl), or by studying other production processes at 
LC 1^, |l^, |l^ and high-energy hadronic colliders^. 

^The most complete list of CP violating dim=6 pmely bosonic operators is given in Concerning them 

we note that the TGC couplings generated by Obw, are identical to those induced by the operators 20B<s>/g or 
20w<s>/g' , defined as the CP violating analogs of Os* and Ow^s respectively. 

^This is particularly needed for the gluon involving operators; (see |l4[). 
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The New Physics contribution from the above operators is described by the effective La- 
grangian 

L-TGC = 9 — —^B<s>+g — —^w<s> + 9 — 2~^w 
l^W 

gg' oiBW , 



where the relations between a^y*) as*, avF; olbw ^ olw ^ and the deviation parameters of Eq.@ 
are given by 

5z = aw*/ (sin 0^ cos x^ = aB$ + api/$ = avi/ 

= — tan^^tiX^ = \^ 

= Q-BW A'y = Q-W 

Kz = — tan^ OwR'y A^ = A^ . (7) 

As it is seen from Eq.(0), the restriction to New Physics generated by dim=6 gauge invariant 
operators, imphes that there are only five independent non-standard triple gauge couplings, 
three of which are CP-conserving |l^] and two CP-violating. 

In order to study the effect of TGC, one usually considers the reaction e+e" — > W^W~, 
taking into account the subsequent decay of the two VF's to a four-fermion final state Q. Such 
final states can be classified in three categories, namely the 'leptonic' ^'i-^^i2^l2, the 'semilep- 
tonic' i'l'e qq' and the 'hadronic' channel ^i^'i^2^2' (where q and q' refer to up- and down-type 
quarks respectively). Semileptonic seems to be the most favoured channel Q for studying TGC, 
since it contains the maximum kinematical information; taking into account that charge-flavour 
identification in a four jet channel is rather inefficient and that the cross section for the leptonic 
mode is suppressed. Thus, in the present paper we study at LC energies, the TGC effect induced 
by the interaction Eq.(|6|) in the processes 

e^e^ i^vi qq (8) 

where ^ is an electron or a muon. The final state TVr qq' is not considered here, due to the 
special difficulties to identify r's in this environment. 

Quite often, the four-fermion final state processes, Eq.(^, are calculated by just taking into 
account contributions from the e^e~ W~^W~ subprocess, which is equivalent to a narrow 
width approximation — > 0. In the classification of the four-fermion production diagrams of 



reference |16|, these graphs are termed as the double-resonant graphs CC03. Such a narrow width 
approximation neglects contributions from single-resonant graphs, which become increasingly 
important at higher energies, (at least in some parts of the phase space) |^]. The 

situation is particularly severe for final states involving , where graphs like the one presented 
in Fig.|^, which involves a t-channel photon exchange, dominate in certain parts of the phase- 
space at higher energies. Moreover, the graph of FigJT] receives contributions from the trilinear 
gauge boson interactions which are not included in the e"*"e~ W~^W~ calculation. In order 
to perform an analysis, as complete as possible, it is therefore mandatory to include in the 
calculation of the processes Eq.(^ all tree-order diagrams, resonant as well as non-resonant 
ones. 
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Figure 1: Single-resonant graph where TGC are contributing. 

Nowadays this is possible, since there exist widely available four-fermion codes, where the 
TGC effects are included beyond the narrow width approximation |^, 17, 18]. In the calculations 
presented in this paper we have used for this purpose, the ERATO Monte-Carlo event generator 
described in [jl^, |2^. The basic ingredients of this generator are the following: 

1. Exact tree-order helicity amplitudes for the processes e^e~ — > ^~^>e, qq' , including all 
trilinear gauge interactions described by Eq.(|l]) |l^, 21 1. 



2. Phase-space generation algorithm based on a multi-channel Monte Carlo approach, includ- 
ing weight optimization [p^ j. 

3. Treatment of the unstable particles contribution consistent with gauge- invariance and high- 
energy unitarity |0, |2^, 24 1 . 



4. Initial state radiation (ISR), based on the structure function approach including soft- 
photon exponentiation as well as hard collinear photon emission in the leading logarithmic 
(LL) approximation up to order O(a^). 

5. Coulomb correction^ to the double resonant (CC03) graphs. 

6. Beamstrahlung effects have also been included via the ^KLgnr]^ algorithm |2 



Apart from the beamstrahlung effects just mentioned, the treatment of the higher order 
corrections in the present study is the same as in the LEP2 case. 

In order to avoid matrix element singularities and to be as close as possible to the experi- 
mental situation, we have applied the cuts 

175° > {01 , Ojet) > 5°, Ee > 10 GeV , Ejet > 10 GeV and mg^g> > 15 GeV . (9) 

Finally, we use the Standard Model input parameters 

Mw = 80.23 GeV, Tw = 2.033 GeV, Mz = 91.188 GeV, Pz = 2.4974 GeV, 

sin2 6'^„ = 0.23103 and a{Mz) = 1/128.07 , (10) 



''For a detailed description see reference 
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while in the ISR structure function a = 1/137.036 is of course used. For the analysis of the 
beamstrahlung effects we have used the TESLA design. 

In order to determine the sensitivity of a given reaction on the TGC parameters one has to 
maximize the likelihood function |28|, whose logarithm is given by 



N 



In Cml = X! lnp(17i, a) 



(11) 



where the sum is running over the event sample under investigation, rjj represents the collection 
of the independent kinematical variables describing the i-th event, the vector a is defined in the 
coupling space as a = {aw^, c^B^ 7 Oiw : CiBW j c^w) i and 



\n=ni 



(12) 



(13) 



is the probability to find an event at the phase-space point Oj. Since the interaction Lagrangian 
is linear with respect to the TGC parameters, one can write the differential cross section in the 
form 

rirr 

(14) 



and similarly the total cross section as 



(15) 



where hatted c's are integrals of unhatted ones over the phase space. 

The sensitivity on the TGC parameters is determined by the so-called information matrix 
given by the second derivative of the likelihood function. 
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(16) 



where 



N 

E[A] = I ll{d^}iPo{ni)}A{n,,...,nN) 

i=l 



(17) 



represents the mean value of a function A{Qi, Oat). If we assume that the maximum of the 
likelihood function is located at a = 0, which reflects the physical expectation that the 'data' 
will be consistent with the Standard Model , the information matrix is given to the lowest order, 
by lij = NBij, with 



^ Co Co / \ Co / \ Co / 



Cl,i 



and 



{A), = J dnpo{n)A{n) , 



(18) 
(19) 
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a ail 



(20) 



a=0 



In the optimal observables approach [30|, equivalent results are obtained by defining the 
phase-space functions 

O, = Sifl (21) 

Coin) 

called optimal observables, whose mean values and covariance matrix determine the sensitivity 
on the TGC parameters. More specifically one writes 

j 

as an unbiased estimator of the components of a, while the corresponding covariance matrix is 
given by 



Via) = ^ • ViO) ■ 



ViO),, = {OiO,)-{0,){0,) , 



{A) = / dn pin) Ain) . 



with ViO) defined by 
and 



Under the assumption that the 'data' are accounted for by the Standard Model, we have that 
ViO) = B, which shows that to the lowest order, the likelihood approach and the optimal 
observables are indeed equivalent. 

Up to now we have considered the so called Maximum Likelihood method. One can improve 
the analysis by considering also the so called Extended Maximum Likelihood (EML). In this 
case we take into account the variation of the total number of expected events as a function of 
the unknown parameters and define 

N 

J~-EML = PNY[pi^i,a) , (22) 



where 

/ AT\N 

-{N) 



and {N) is the number of expected events. The analysis proceeds as before and the correlation 
matrix is expressed as 

= /£M£Li\ (23) 

\ Co Co /o 

while the information matrix becomes lij = {N)Bij. The corresponding optimal observables in 
the EML approach are 

0. = n'-^. (24) 

Coin) 

In the sequel we perform one- as well as five-dimensional investigations. One-dimensional 
investigations assume that all but one of the a^'s are non- vanishing, and lead to parameter errors 
(Isd) given by 
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The multidimensional case where all five NP couplings are considered simultaneously, is also 
treated. Diagonalizing then the symmetric correlation matrix B, one first determines its eigen- 
values Xi and eigenvectors , normalized so that • ej = Sij . For each Cj , the parameter 

af = Ci - a (26) 

is then defined, for which the (Isd) error is given by 

In all calculations presented in this paper, N is taken to be the predicted Standard Model 
number of events defined by A'' = 4 L (Tsm-, where asM is the corresponding total cross section, 
L is the integrated luminosity and the factor 4 takes into account the four equivalent channels 
described by the same matrix elements; i.e. e^e~ — ^ (-"ve. ud, t'vii cs, i'^i'e du and i'^i'e sc. 



L = 50 fb"' 
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Figure 2: The one (solid) and two (dashed) standard deviations limits on aB^-ctw^, using the 
optimal observables and the EML methods for -s/s = 800 GeV and unpolarized beams. In the 
upper part optimal observables and EML are hardly distinguishable, whereas in the lower part 
EML exhibits a secondary minimum. 



At this point we address the question, how accurately the optimal observables approximation 
describes the ML (or EML) function. To this end we calculate the In Cm L by replacing the sum 
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appearing in Eq.(ll) by an integral over the expected probability distribution which is assumed 
to be the one predicted by the Standard Model 



N 

In Cm L = N j JJl^^^^i Po(^^i)} a) • (28) 

i=l 

It is clear that optimal observables and ML methods become identical in the limit — > oo, 
since then only the leading term in the expansion of the likelihood function survives; and this 
is exactly the term proportional to the information matrix. On the other hand for relatively 
low statistics, the nonlinearity of the likelihood function becomes important and the optimal 
observables approximation breaks down. These features are shown in Fig.|2|, where the one, 
and two standard deviation limits on (a_B4>) oiw<^) are considered, for the muon channel and 
y/s = 800 GeV. In the upper part of the figure the value of the integrated luminosity is taken 
to be -L = 50 fb~^, which is the expected nominal value, whereas in the lower part a much lower 
luminosity, L = 5 fb~^, has been used. 

We have checked that for all nominal LC energies and luminosities, the optimal observables 
approximation gives identical results to those of the conventional likelihood approach. Fur- 
thermore the fact that the expected sensitivities on the TGC parameters are predominantly 



determined by the linear terms in the expansion of the differential cross section, Eq.(14), shows 
the self-consistency of our original assumption that a parameterization of TGC in terms of dim=6 
operators should be adequate. 

On the other hand, from the point of view of a weighted Monte-Carlo approach, which is 
frequently used in the phenomenological analyses, the optimal observables method offers a very 
efficient fast and economic way to estimate not only the sensitivity of a given process on a single 
TGC parameter (or any kind of 'deviation' parameter), but also their full covariance matrix, 
which is of great importance for multiparameter analyses. 

Finally we would like to mention that the general experimental problem of how to overcome 
the ISR as well as the detector resolution induced difficulties in the reconstruction of the event 
kinematics, is not addressed here. We only note that this problem exists also for the current 
LEP2 experiments and that detailed experimental simulations at linear collider energies can be 



found in the proceedings of the DESY-ECFA Workshop on Linear Colliders |27|. Moreover our 
experience from LEP2 studies Q shows that, despite the abovementioned problem, a very good 
estimate of the sensitivity on the TGC can be obtained by an analysis of the kind used in our 
present study. 

In Table |^ we present the results for the correlation matrix Bij involving all CP-conserving 
and CP-violating couplings, at 500 GeV center of mass energy. The total cross sections are 
also presented. As is evident from this table, the correlations between the different Oj's are not 
negligible in general, which suggests that an analysis taking them into account, is indispensable. 

An other very interesting result, is that electron and muon channels exhibit a complementary 
behaviour: electron channel gives the highest production rate, which means a better statistics, 
whereas the muon-channel exhibits a higher sensitivity on TGC. 

In Table ^ we show the eigenvalues of the correlation matrix, as well as the correspond- 
ing eigenvectors. These eigenvectors define directions in the five-parameter space, which are 
uncorrelated, so that parameter errors can be safely extracted. This table shows that for the 



8 





Vlr + Vbl 


Vrl 


e 


14.53 


4.87 


3.26 - 


-0.0032 


-0.27 


9.66 


-10.20 


0.33 


0.59 


0.027 




4.87 


3.67 


0.33 


0.016 


-0.072 


-10.20 


70.33 


-0.11 


-1.80 


-0.096 




3.26 


0.33 


20.75 


-0.011 


0.38 


0.33 


-0.11 


3.50 


0.012 


-0.016 




-0.0032 


0.016 


-0.01 


0.28 


-0.78 


0.59 


-1.80 


0.012 


6.66 


0.22 




-0.27 


-0.072 


0.38 


-0.78 


21.53 


0.027 


-0.096 


-0.016 


0.22 


3.46 




456(5) 


207(5) 




23.33 


7.33 


5.66 


0.013 


-0.060 


258.42 


-513.87 


12.94 


-1.35 


-0.12 




7.33 


5.25 


0.71 


0.033 


-0.054 


-513.87 


2478.22 


-9.50 


12.36 


0.28 




5.66 


0.71 


33.00 


-0.032 


0.29 


12.94 


-9.50 


2.71 


0.27 


-0.0013 




0.013 


0.033 


-0.032 


0.28 


-1.33 


-1.35 


12.36 


0.27 


178.16 


12.11 




-0.060 


-0.054 


0.29 


-1.33 


33.85 


-0.12 


0.28 


-0.0013 


12.11 


2.67 




267(2) 


6.01(4) 



Table 1: The correlation matrix for e and fi channels at 500 GeV for the TGC parameters 
aw^, ob^, ct-w ^olbw lOiw ■ Also shown are the cross sections as well as their Monte Carlo errors 
in femptobarns. Here VrlCPlr) corresponds to e^e^ (^Z^ij) initial state polarization. 



unpolarized beams case, the dominant eigenvalues correspond to directions in the five-parameter 
space related predominantly to aw, Oiw^, and aw, whereas the lowest ones are related to as^ 
and asw- The picture becomes almost opposite in the case that the electron (positron) beams 
are right- (left)-polarized. 

As far as the polarization is concerned, we see that passing from unpolarized to right- 
polarized electron-beam, results to a much higher sensitivity for as$ and obw couplings. These 
phenomena are much more pronounced for the muon-channel. This effect, which has been also 
observed in on shell studies of e~e^ — > W~W^ [^], reflects the fact that different TGC param- 
eters contribute to different helicity amplitudes, especially in the high-energy regime. In Fig.^ 
we show how the information from both polarizations can, in principle, be used to disentangle 
different TGC parameters, based on the fact that the corresponding covariance matrices are 
very different. 

Finally in Table |, one-standard-deviation errors are presented by combining electrons and 
muons as 



A^ = 4L(cj(^)+a('^)) . 

In our studies, L is taken to be 20 fb'^ at 500 GeV, 10 fb^^ at 360 GeV and 50 fb^^ at 800 
GeV. For the results concerning polarized beam scattering, we used 

^polarized ~ 4 -^unpolarized " 
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36.10 -0.460 -0.129 -0.872 0.004 -0.099 
33.90 0.003 0.002 -0.007 0.039 -0.999 
22.88 0.816 0.320 -0.479 0.0005 0.018 
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Vrl 


72.06 0.161 -0.986 0.002 0.028 0.001 
8.07 -0.963 -0.163 -0.067 -0.198 -0.011 
6.56 -0.199 -0.0043 -0.018 0.977 0.068 
3.48 0.062 0.007 -0.894 -0.034 0.440 
3.44 -0.028 -0.003 0.440 -0.059 0.895 


2591.53 0.215 -0.976 0.004 -0.005 -0.0001 
178.98 -0.040 -0.003 -0.004 -0.996 -0.068 
145.96 0.973 0.214 0.073 -0.040 -0.003 
1.87 0.068 0.010 -0.935 -0.022 0.346 
1.83 -0.024 -0.003 0.346 -0.064 0.935 



Table 2: The eigenvalues (2nd and 8th columns) and the corresponding eigenvectors of the 
correlation matrices given in Table 



In order to study the effect of the correlations among the TGC parameters we distinguish two 
cases: 

• The Id-case is based on the very strong and often made assumption that only one of the 
TGC parameters {aw^, as*, Oiw-, o^bWi o:w) is non vanishing. A very contrived form of 
NP at high energy scale is needed, in order to create such a situation where only one of the 



operators appearing in Eqs.(^,§) is generated at low energies |^. This case corresponds 
to the 'one-dimensional log-likelihood fit'. 

In the 5(i-case, on the contrary, the errors are calculated according to Eq.(p7|), where the 
full correlation matrix is taken into account and no a priori assumption on the size of 
the parameters has been made. Although in this case the presented errors correspond to 
directions in the five-parameter space defined by Eq.(p6D, which are not generally identical 
to the ones defined by the original parameters, we kept the same notation, since the former 
are rather close to the latter: for instance a^,^ is mainly composed by aw^ and so on for 
the other TGC parameters p3|. 



As it is seen from Table |3|, moving from the Id-case to the 5d-case, the change on the one- 
standard-deviation errors reach the level of 40%. Moreover the less sensitive the TGC parameter 
is, the more the correlations affect its error. It should be mentioned however that the correla- 
tions among the different TGC parameters do not dramatically change the order-of-magnitude 
estimate of their sensitivity based on single-parameter considerations. 



Concerning the CP violating interactions, we should note that in reference |34], bounds on 
the CP-violating couplings k-^ and A.y have been derived, on the basis of their contribution to 
the electric dipole moment (EDM) of the meutron. Besides the fact that these bounds depend 
on several details, the most they imply is a strong linear relation between and X-y. Direct 
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Figure 3: The one (solid) and two (dashed) standard deviations limits on as#-aiy$ for unpolar- 
ized and polarized (e^ej) beams. The inner lines are from electron channel whereas the outer 
ones are from muon channel. 



measurements of these couplings, as well as of their Z-boson counterparts, will therefore be useful 
because they will provide detail information on the whole CP-violating TGC parameter-space. 

In studying the sensitivity on the TGC one usually neglects possible correlations with other 
electroweak parameters like for instance the non-standard contributions to V f f vertices, where 
V stands for Z or W . This is rather well justified because the latter are usually much more 
constrained than the former as it is indeed the case at LEP2, where TGC determination is 
expected to reach the level of 0.01 to 0.1, compared with the constraints on the V f f couplings 
coming mainly from LEPl analysis, which are at the level of 10~^ [^]. On the other hand, as 
it is also suggested by our analysis, at LC energies the TGC can be tested to a precision much 
higher than that of LEP2, and therefore it becomes interesting to study the correlations among 
the TGC and the other electroweak couplings 



1,0 



Finally we would like to mention that a detailed comparison of our study with those presented 
in references |^, 32] based on the on-shell production e~e^ — > W~W^ is not possible, due mainly 



to the fact that we are using different analysis methods. Nevertheless both approaches agree in 
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Table 3: One standard deviation errors on TGC parameters. At 500 GeV we show also the 
effect of the correlations between different five TGC parameters in the bd case for unpolarized 
beams, whereas in parentheses we show the corresponding errors from a right-handed polarized 
(left-handed) electron (positron) beam, as explained in the text. 

the order-of-magnitude estimate of the expected sensitivity on the TGC. 
We conclude by summarizing the results of our study: 

1. We have presented a five-parameter description of the non-standard trilinear gauge cou- 
plings, which includes all gauge-invariant contributions to the lowest order. We then 
analysed their contribution to the scmileptonic reactions e~^e~ — ^ ^~i^e qq' for £ = e and 

i = fi at LC energies, and showed that a measurement of all five parameters is possible, 
with a sensitivity covering a rather wide range starting from 1.1 x 10~^ (Isd) at 360 GeV 
for Q!_Biy, (worst case), and going down to 3 x 10"'^ (Isd) at 800 GeV for aw and aw, 
(best case). 

2. The electron channel, due to the onset of the singlc-VF production mode, gives the dom- 
inant contribution to the total cross section at LC energies, whereas the muon channel 
exhibits a higher sensitivity on the TGC parameters. Therefore their overall contribution 
to the error on the TGC determination become equally important. 

3. Polarization effects are important in order to disentangle different TGC contributions, 
leading to a substantial improvement of the sensitivity on the TGC parameters, especially 
for aB^ and olbw- 
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